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Donors in silicon hold considerable promise for emerging quantum technologies, due to the their 
uniquely long electron spin coherence times. Bismuth donors in silicon differ from more widely stud- 
ied Group V donors, such as phosphorous, in several significant respects: they have the strongest 
binding energy (70.98 meV), a large nuclear spin (/ = 9/2) and strong hyperfine coupling constant 
[A — 1475.4 MHz). These larger energy scales allow us to perform a detailed test of theoretical 
models describing the spectral diffusion mechanism that is known to govern the electron spin deco- 
herence of P-donors in natural silicon. We report the electron nuclear double resonance spectra of 
the Bi donor, across the range 200 MHz to 1.4 GHz, and confirm that coherence transfer is possible 
between electron and nuclear spin degrees of freedom at these higher frequencies. 



Electron and nuclear spin coherence of donors in sili- 
con is of great importance for a number of proposals for 
Si-based quantum technologies |IH1] ■ These schemes cite 
among the advantages for Si donor quantum bits: long 
coherence times (exceeding tens of milliseconds for the 
electron and seconds for the nucleus in the case of Si:P), 
high-fidelity manipulation through a combination of mi- 
crowave and radiofrequency pulses and integration 
within silicon devices for measurement [3HS]- Magnetic 
resonance studies on P-donors in Si have examined elec- 
tron spin coherence in natural Si [lOl lll| and its depen- 
dence on increasing ^^Si concentration in the host [T^ . 
as well as the storage of coherent electron spin states in 
the ^^P nuclear spin [13j . 

Although there has been a focus on the P-donor in Si, 
other Group V donors such as Bi also possess attractive 
qualities as quantum bits [Hj. There have been rela- 
tively few recent studies on the Si:Bi, an exception be- 
ing a photoluminescence study showing dynamic nuclear 
polarization of the ^"^Bi through optical pumping [T5]. 
Bismuth is the deepest group V donor, with a binding 
energy of 70 meV [TB] and the largest Group V hyperfine 
coupling of 1.4754 GHz to the 1 = 9/2 nuclear spin of 
209gj These parameters differ substantially from the P 
donor (44 meV and 117.52 MHz), raising the question of 
whether the same decoherence mechanisms and method- 
ology for nuclear spin manipulation are applicable. In 
this letter we examine spin decoherence of Bi donors in 
natural silicon, as well as electron-nuclear double reso- 
nance (ENDOR) spectroscopy to probe the transitions of 
the ^°^Bi nuclear spin. We find that these measurements 
compare well with P donors in natural silicon, providing 
strong motivation to pursue ^'*Si:Bi material. 

Natural Si:Bi samples were obtained from from ul- 
trapure "'^^Si starting material by a floating-zone tech- 
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FIG. 1: (Colour online) Electron spin resonance of Si:Bi 
donors at X-band (9.7 GHz), a) Energy levels of the coupled 
electron-nuclear spin system as a function of magnetic field. 
The allowed EPR transitions at 9.7 GHz are indicated with 
vertical lines, each corresponding to a different mj projection 
for the 7 = 9/2 ^°^Bi nuclear spin, b) Experimental elec- 
tron spin echo intensity as a function of magnetic field yields 
the 10 expected resonances. Temperature=16 K, simulation 
parameters: A = 1475.4 MHz, = 2.00, 7Bi =6.962 MHz/T. 



nique, as described in Ref |18j and had a room tem- 
perature resistivity of 4.5 llcm implying a Bi concen- 
tration of 10^^ cm~^. Pulsed EPR measurements were 
performed using a Bruker Elexsys 680 X-band spectrom- 
eter, equipped with a low temperature helium-fiow cryo- 
stat (Oxford CF935). A TWT amplifier was used to 
provide EPR tt/2 and tt pulses of 20 and 40 ns respec- 
tively. For ENDOR measurements an Amplifier Research 
60 W solid state CW amplifier (0.8-4.2 GHz) and an ENI 
100 W (1.5-400 MHz) were used, depending on the fre- 
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FIG. 2: (Colour online) Two-pulse electron spin echo decay of 
the Si:Bi donor as a function of angle of the applied magnetic 
field Bq with respect to the [100] crystal axis. Crystal rotation 
is performed in the [100]-[011] plane. T = 12 K, Bo = 5663 G. 

quency range, driven by an Agilent PSG Analogue Sig- 
nal Generator. RF pulses in the range 4-15 fis were used 
for a TT-pulse on the ^°^Bi nuclear spin, depending on 
frequency. EPR/ENDOR spectra were simulated using 
Easyspin p^ . 

The Si:Bi electron/nuclear spin system can be de- 
scribed by an isotropic spin Hamiltonian (in angular fre- 
quency units): 

no^L^eS,-uiI,+A-S-I, (1) 

where uJe = gl^Bo/h and w/ = gi/3nBo/h are the elec- 
tron and nuclear Zeeman frequencies, g and gi are the 
electron and nuclear g-factors, /3 and (3n are the Bohr 
and nuclear magnetons, h is Planck's constant and Bq is 
the magnetic field applied along z-axis in the laboratory 
frame. The donor electron spin S=l/2 (g — 2.0003) is 
coupled to the nuclear spin / = 9/2 of ^°^Bi through a 
hyperfine coupling A = 1475.4 MHz |I7]. At high mag- 
netic fields (i.e. We ^ A{I + 1/2)), this leads to ten 
equally spaced resonances in the EPR spectrum, each 
corresponding to a transition Ams = ±1 for a given mj 
projection. As shown in Figurejl] measurements made at 
X-band (9.7 GHz) are not entirely in this high-field limit. 
The EPR spectrum of Si:Bi was recorded by monitoring 
the electron spin echo (ESE) intensity as a function of 
magnetic field. A linewidth of 4 G was measured for 
each of the ten EPR lines, consistent with inhomoge- 
neous broadening from unresolved hyperfine coupling to 
the surrounding ^^Si nuclear spins (~ 5% natural abun- 
dance). 

Electron spin decoherence of P-donors in natural sil- 



icon is known to be dominated by spectral diffusion, a 
mechanism in which spin flip-flop of surrounding ^^Si 
nuclei modulates the electron Zeeman energy through 
both contact and dipolar hyperfine coupling [501 HI]- The 
^^Si nuclei closest to the donor are a 'frozen core' which 
do not flip-flop due to the strong spatial dependence of 
their coupling to the donor electron spin, causing detun- 
ing between adjacent ^^Si nuclear spins and suppressing 
nuclear flip-flop transitions that are allowed in the bulk 
material. Similarly, those furtherest away are too weak 
to influence the donor electron. There is therefore an 
'active shell' in the region where the dipolar coupling 
between neighbouring ^^Si spins is comparable to their 
coupling to the donor electron spin, which is responsi- 
ble for spectral spectral diffusion. This mechanism has 
been predicted [221 US] and shown experimentally [TT] to 
have an angular dependence corresponding to the depen- 
dence of the dipolar coupling between nearest-neighbour 
^^Si spins on the crystal orientation with respect to the 
applied magnetic field. Figure [2] shows the measured elec- 
tron spin echo decay traces as a function of angle, mea- 
sured at the high-field line (5663 G, to/ = -9/2). The 
behaviour is qualitatively similar to that of P-donors in 
""^'Si, though decay times are approximately 30% longer. 
This can be primarily attributed to the greater binding 
energy of the Bi donor compared with P, which shrinks 
the effective Bohr radius of the Si:Bi donor in comparison 
to Si:P, reducing the size of the 'active region'. There is 
also a secondary effect resulting from the much stronger 
hyperfine coupling to the donor nuclear spin, as described 
below. 

Although other Group V donors are well into the high- 
field approximation at X-band (9.7 GHz, 0.35 T), the 
large hyperfine coupling to ^"^Bi causes some level mix- 
ing, as shown in Figure [T] As a measure of the sensitiv- 
ity of a transition frequency / to a change in magnetic 
field B, we can extract an effective gyromagnetic ratio 
7eff = df/dB which differs substantially from that of a 
free electron (7 = 28.0 GHz/T), and also varies for each 
of the hyperfine lines across the EPR spectrum: from 
7eff - 26.7 GHz/T (for to/ = ±9/2) to 18.4 GHz/T 
(for TO/ — ±1/2). The resulting change in alters the 
coupling of the donor electron to surrounding ^^Si and 
thus selects a different 'active shell', with different statis- 
tics of pairwise ^^Si coupling. Thus, though the dipolar 
coupling between neighbouring ^^Si spins is primarily re- 
sponsible for setting the timescale of spectral diffusion, 
we may expect an effect due to varying 7off . Through sim- 
ulations based on cluster expansion technique of Ref [55] 
we have calculated the effect of varying jeS, which pre- 
dict a ^ 5% increase in spectral diffusion times, TgD, as 
measured on the to/ = ±1/2 hyperfine line, compared 
to TO/ = ±9/2. The simulations use a Kohn-Luttinger 
wavefunction with the 70 meV binding energy for Bi. 

We fit the echo decay traces, such as those in Figure [2] 
to a combination of an orientation-independent T2, com- 
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FIG. 3: (Colour online) Extracted spectral diffusion times, 
TsD, as a function of crystal orientation and magnetic field, a) 
A fit to the electron spin echo decay curves provides a measure 
of Tso as function of angle of the applied magnetic field Bq 
with respect to the [100] crystal axis, performed on two of the 
ten hyperfine lines: 5663 G (blue, circles) and 2542 G (red, 
squares). The solid curve is a spline fit through the data for 
5663 G, which is multiplied by 1.05 to produce the dashed 
curve, b) Inset shows the magnetic field dependence of Tsd 
for six field positions (6 = 0°), showing good agreement with 
the results of the simulation (see text). Temperature — 12 K. 



bined with an orientation-dependent Tsd , through an ex- 
pression of the form: 



(2) 



The spectral diffusion mechanism has a characteristic 
stretched exponential coherence decay, with a typical 
value of n between 1-4 depending on the regime of spec- 
tral diffusion [531 I^^^ES] . The stretching factor n was 
found to be independent of crystal orientation, while Tsd , 
plotted in Figure[3]shows the expected orientation depen- 
dence with a maximum when the applied magnetic field 
Bq is oriented along [100], and a minimum when oriented 
along [111]. Comparing the values measured at two hy- 
perfine lines: mi = -1/2 (2542 G) and m/ = -9/2 
(5663 G) we see agreement with the j^s dependence 
predicted by the simulations. The stretching factor n 
showed no significant field dependence: 7125420 = 2.30(7) 
and n5663G — 2.34(5), consistent with the simulations 
which predict 2.30(1). 

The inset of Figure |3] shows the predicted values of 
Tsd from the simulations, which are within ^ 20% of 
the experimental values. The magnetic field dependence 
of Tsd is also well represented, showing a maximum for 
TO/ = ±1/2. The simulations were performed using the 
cluster correlation expansion |27j but it is sufficient to use 
a simple pair approximation that includes effects from 
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(Colour online) Electron nuclear double resonance 
^"'■'Bi in silicon, a) Dashed curves show the the- 
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oretical ENDOR frequencies as a function of field, for each 
Ami = 1 transition. Symbols represent those values mea- 
sured by Davies ENDOR at X-band, at each of the 10 reso- 
nant fields of the EPR spectrum, b) and c) show two typical 
ENDOR peaks, at 1345.4 and 677.4 MHz with RF yr-pulse 
length = 4 and 7 jj,s, respectively. 



each pair of nuclear spins independently |25l I28j . 

The fitting error in the residual T2 parameter is large 
when it is much longer than Tsd, however, our extracted 
values of a few milliseconds are consistent with being lim- 
ited by electron spin relaxation (Ti) at this temperature 
(12 K). We find that Ti is well described by a first-order 
Raman mechanism (T^'') in the temperature range 8- 
16 K, as proposed by Castner for the range 19-25 K [29] 
(see Supporting Information). 

We now turn to an investigation of the ^°^Bi nu- 
clear spin transition through electron-nuclear double res- 
onance (ENDOR). We used the Davies ENDOR sequence 
("■mw — ■"'rf — 7r/2,ji„ — r — TTya^ — T — ccho) at each of the ten 
hyperfine lines to map out the set of 36 distinct ENDOR 
lines observable at a given microwave frequency: for each 
hyperfine line (i.e. EPR transition) there are four observ- 
able ENDOR transitions (Ato/ = 1), apart from at the 
ends of the EPR spectrum (to/ = ±9/2) where there are 
only two. The measured frequencies are shown in Fig- 
ure [4^, along with theoretical curves showing ENDOR 
frequencies as a function of magnetic field. 

Two typical ^"^Bi ENDOR lines are shown in Fig- 
ure |4)d and|4j;. The ENDOR hnewidth ranges from 100- 
700 kHz, depending on the particular transition and field 
being measured. In general, the linewidths are broader at 
lower magnetic fields or rf frequencies, except for the two 
lines corresponding to the to/ — —9/2:— 7/2 transition 
which are about 350 kHz wide. The ENDOR linewidths 
we observe (plotted in the Supporting Information) are 
well described by a combination of two factors: At lower 
magnetic fields (< 400 mT), the linewidth arises from 
the random dipolar field of surrounding ^^Si, as in the 
case for the EPR linewidth. This effect is directly re- 
lated to the gradient of the field/frequency curves shown 
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in Figure [4^ which flatten out as the high-field approxi- 
mation becomes valid (this also accounts for the greater 
linewidth of the m/ = —9/2:— 7/2 transition). Corre- 
spondingly, this broadening mechanism is not significant 
for ENDOR in other Group V donors at X-band given 
their much weaker hyperfine couplings. Instead, ENDOR 
linewidth in such donors arises from an inhomogeneity in 
the hyperfine coupling to the donor nucleus, due to a 
variation in the dielectric constant of the material within 
the donor wavefunction caused by the random distribu- 
tion of ^^Si. Such a mechanism could be responsible for 
the ENDOR linewidth in Si:Bi at higher magnetic field 
(< 400 niT) and from our measurements we can put 
an upper bound of the inhomogeneity in a(^''^Bi) to be 
(< 0.02%). 

It has been shown that the ^^P donor nuclear spin can 
provide a valuable resource for storing the coherent state 
of the electron spin for times exceeding seconds [13] . The 
larger nuclear spin (/ = 9/2) of ^'^^Bi provides a corre- 
spondingly larger Hilbert space for storing electron spin 
qubits, though also introduces more potential relaxation 
mechanisms for the nuclear spin. As we have found for 
Si:P, high-fidelity storage/retrieval of the electron coher- 
ence requires narrow EPR and ENDOR lines and thus 
a Si-enriched host material. We have investigated the 
potential of the •^'^^Bi nuclear spin for quantum mem- 
ory, and found that as expected, the '^'^^Si host limits the 
store/retrieve fidelity to ~63% (see Supporting Informa- 
tion) . We can nevertheless measure the decay of nuclear 
coherence and found that it is limited by the effect of 
Tic (first order Raman) processes at temperatures above 
10 K — random fiuctuations in the electron polarisation 
drive decoherence of the strongly coupled nuclear spin. 
Below this temperature, is limited to about 15 ms. 
Further work using ^^Si:Bi will be required to explore this 
limit and investigate the sources of nuclear decoherence. 

We have found electron spin decoherence of Bi donors 
in natural Si to be dominated by the same spectral diffu- 
sion mechanism found in the case of Si:P, however thanks 
to the smaller Bohr radius of Bi the effect is weaker than 
for P, leading to 30% longer T2 times. Despite the high 
ENDOR frequencies necessary to probe the ^°^Bi nuclear 
transitions, it is possible to excite each transition with a 
fidelity determined by the ENDOR linewidth of a few 
100 kHz allowing us to demonstrate the possibility of 
storage and retrieval of electron spin coherence in the 
^°^Bi nuclear spin. Applying such techniques to ^^Si:Bi, 
we would anticipate the ability to store and retrieve mul- 
tiple electron spin qubits with high fidelity within the nu- 
clear spin. Finally, we note that the large energy splitting 
present at zero applied magnetic field makes the Si:Bi 
donor spin an attractive candidate for coupling to super- 
conducting resonators. 

Note added: Interesting experimental and theoretical 
investigations of Si:Bi in the context of quantum infor- 
mation are reported in parallel studies [301 131] • 
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SUPPLEMENTARY INFORMATION 
Spin- lattice relaxation times (Tic) 

Electron spin-lattice relaxation times Tic were 
recorded as a function of temperature between 8 and 
25 K using a three-pulse inversion recovery sequence 
(tt)^ — t — (7r/2)_|_^ — (tt)^ — (echo) and two-step phase 
cycling. The inversion recovery plots were well fit to a 
mono-exponential decay and temperature was recorded 
with a calibrated thermometer below the microwave res- 
onator. 



SI 



Tie is shown as a function of temperature in Fig 
together with the least-squares fit to a model of the form 
Tie = AT-^ where A = 4.3 (±0.3) x 10'* sK^, consistent 
with the relaxation being driven by an inelastic Raman 
phonon scattering process [Tj . The Tie data was collected 
from the mj — —9/2 line at 577.2 mT; we found Tie to 
be independent of the choice of mi . 



Si:Bi ENDOR Linewidths 

Linewidths of the transitions between nuclear sub- 
levels were recorded using a Davies ENDOR sequence: 



TTrf - 7r/2n 



T — echo. The RF pulse 



power was optimised by recording nuclear Rabi oscilla- 
tions while the RF pulse length was chosen to avoid in- 
strumental broadening of the ENDOR lines. ENDOR 
transition frequencies and linewidths were recorded on 
each of the ten EPR lines and the resulting distribution 
of linewidths is displayed in Fig. |S2[ together with the- 
oretical calculations of the expected line broadening due 
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FIG. SI: Temperature dependence of spin-lattice relaxation 
time in Si:Bi between 8 and 25K, showing a Raman-like T"'^ 
behaviour. 



to (a) strain in the hyperfine coupling to ^°^Bi and (b) 
random magnetic fields due to dipolar coupling to ^^Si 
nuclear spins. The observed increase in linewidth at lower 
magnetic fields is consistent with the latter mechanism, 
though A-strain may contribute to ENDOR linewidth at 
fields above about 400 mT. 



Fidelity of Si:Bi Nuclear Memory 

Using the method described in Ref , we demonstrate 
the ability to store some of the electron spin coherence in 
the nuclear spin, and subsequently recover it, as shown 
in Fig. [83] It is not possible to fully excite the ENDOR 
hyperfine line, leading to a much lower transfer fidelity 
than observed for P-donors in isotopically purified ^*Si, 
but comparable to that observed for "^'Si:P. 

Phase cycling of both microwave and rf pulses is used 
to confirm the recovered signal is a result of storage in the 
nuclear spin and is not an artefact such as a stimulated 
echo. 

Comparing the intensity of the recovered electron spin 
echo with that which was stored allows for a rough as- 
sessment of the performance of the transfer — in this 
case yielding /stored ~ 25%. 

To compare the transfer fidelity with other quantum 
memory experiments, we estimate the stored and re- 
covered density matrices and compute the fidelity of 
that retrieved with respect to that stored using J- = 
Tr (profPstorcd)- Wc ucglect the az component in both 
cases, and use the integrated electron spin echoes in the 
real and imaginary channels of the quadrature detector 
to extract (Tx and ay components. As in Ref [2, the ref- 
erence echo is designated as a pure state (i.e. no identity 
component) and used to normalise the intensity of the 
recovered echo. 

By this method, we estimate the fidelity of the 
store/retrieve process to be = (1 + /stored) /2 = 0.635. 
To avoid over-interpreting this result, it should be noted 
that recovering zero electron coherence (i.e. retrieving 
only the identity matrix) would yield a fidelity of 0.5. 



* Electronic address: |john.mo rton@materials.ox.ac.uk| 
[1] T. G. Castner. Phys. Rev., 130(l):58-75, Apr 1963. 
[2] J. J. L. Morton, A. M. Tyryshkin, R. M. Brown, 
S. Shankar, B. W. Lovett, A. Ardavan, T. Schenkel, 
E. E. Haller, J. W. Ager, and S. A. Lyon. Nature, 
455(7216) :1085-1088, Oct 2008. 
[3] A. Schweiger and G. Jeschke. Principles of Pulse Elec- 
tron Paramagnetic Resonance. Oxford University Press, 
Oxford, UK ; New York, 2001. 



7 



dcD / dAgi dcD / dB Observed 




200 400 600 200 400 600 200 400 6( 

Magnetic Field (mT) Magnetic Field (mX) Magnetic Field (mX) 



FIG. S2: ENDOR linewidths of hyperfine transitions in the Si:Bi systems observable using a microwave frequency of 9.7 GHz: 
(a) (A) predicted due to a 0.23 MHz strain in the ^"^Bi hyperfine parameter A (corresponding to 0.02%), (b) (o) predicted 
due to random variation in magnetic field of ~ 2 G. (c) (□) Experimental data. 
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FIG. S3: Electron spin echo transients in Si:Bi comparing a reference (Ref) electron spin coherence with one that is stored in the 
nuclear spin and then recalled. The loss in signal intensity between the Reference and Recalled echoes is due to the imperfect 
fidelity of the coherence transfer process. For both cases, the real and imaginary components of the spin echo, measured using 
a quadrature detector, are shown separately. These signals are proportional to measurements of the ax and ay components of 
the density matrix, respectively |3]. Temperature = 14 K 



